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Abstract 

A  simple  method  is  presented  for  estimating  the  combined  influence  of 
strain-hardening  and  strain-rate  sensitivity  on  the  permanent  deformation 
of  rigid-plastic  structures  loaded  dynamically.  A  study  is  made  of  the 
particular  case  of  a  beam  supported  at  the  ends  by  immovable  frictionless 
pins  and  loaded  with  a  uniform  impulse.  The  results  of  this  work  indicates 
that  considering  strain-hardening  alone  when  appropriate  or  strain-rate 
sensitivity  alone  gives  permanent  deformations  which  are  similar  to  those 
predicted  by  an  analysis  retaining  both  effects  simultaneously. 
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strain-rate  sensitivity  coefficient  defined  in 
equation  (38) 

modulus  of  elasticity 

thickness  of  beam 

coefficients  given  by  equations  (64)  and  (65) 
semi-length  of  beam 

longitudinal  bending  moment  per  unit  width  of  beam 

o  H2/4 
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axial  force  per  unit  width  of  beam 

o  H 
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shear  force  per  unit  width  of  beam 
radius  of  curvature  of  mid-plane  of  beam 
time 

initial  velocity  of  beam 
-k  sin$ 

constants  defined  in  equations  (61),  (76) 

uniform  distributed  pressure  per  unit  area  of 
undeformed  beam 

M/M 

o 

N/N 

o 

dimensionless  membrane  forces  defined  by  equations 
(80,  84,  79,  74,  85,  59) 

strain-rate  sensitivity  coefficient  defined  in 
equation  (38) 
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Notation  (continued) 
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ratio  of  the  slopes  of  the  elastic  and  plastic 
portions  of  the  stress-strain  curve 

distance  measured  along  deformed  mid-plane 

time 

duration  of  first  stage  of  motion 
t1/L2 

time  at  which  beam  reaches  its  permanent  position 

displacements  defined  in  Fig.  1 
maximum  permanent  transverse  displacement 

distance  defined  in  Fig.  1  (measured  from  beam  center) 
distance  measured  from  left-hand  support  of  beam 
distance  defined  in  Fig.  1 

2V  H  1/p 
(-2-) 

DL2 


CD 

3  1  T 

k=l,3,5  k1* 

y  vH2/L2 

61 .  Vi 

e  axial  strain  of  mid-plane 

ez  axial  strain  at  distance  z  from  mid-plane 

defined  by  equation  (77) 
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curvature 


Notation  (continued) 
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radius  of  a  traveling  hinge 
stress 

yield  stress  in  simple  tension 
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slope  of  the  mid-plane  of  beam 


k=l,3,5  k3 
fij/H 


difference  between  the  values  of  the  considered 
quantity  on  either  side  of  a  traveling  hinge 


1 .  Introduction 


Parkes  [1]  examined  the  permanent  deformation  of  cantilever  beams  loaded 
dynamically  and  observed  that  a  simple  rigid,  perfectly  plastic  analysis  over¬ 
estimated  considerably  the  final  maximum  deflections.  The  rather  significant 
discrepancies  between  experimental  results  and  theoretical  predictions  were 
accounted  for  by  considering,  in  an  approximate  manner,  the  influence  of 
strain-rate  sensitivity  on  the  dynamic  plastic  bending  moment.  Ting  [2]  ana¬ 
lyzed  a  rigid,  perfectly  plastic  cantilever  beam  loaded  dynamically  at  its  tip 
and  indicated  that  geometry  changes  when  treated  rigorously  were  responsible 
for  part  of  the  discrepancies  between  Parke's  theory  and  experiments  [1]. 

Bodner  and  Symonds  [3,4]  conducted  more  exhaustive  tests  on  cantilever  beams 
loaded  dynamically  and  observed  that  strain-hardening  was  not  very  important 
while  strain-rate  sensitivity  must  be  considered  throughout  the  entire  defor¬ 
mation  history.  Ting  and  Symonds  [5]  analyzed  the  plastic  deformation  of  a 
cantilever  beam  with  an  attached  tip  mass  which  was  subjected  to  a  rapid 
transverse  velocity  change  at  the  base  and  found  that,  if  the  strain-rate  de¬ 
pendence  of  yield  stress  and  geometry  changes  were  considered,  then  the  predic¬ 
tions  showed  good  agreement  with  corresponding  experimental  results. 

Parkes  [6]  developed  the  earlier  theoretical  work  of  Lee  and  Symonds  [7] 
in  order  to  describe  the  behavior  of  a  rigid,  perfectly  plastic  encastr«  beam 
struck  transversely  at  any  point  on  the  span.  The  predictions  were  compared 
with  some  experimental  values  recorded  on  steel,  brass,  and  duralumin  beams, 
the  supports  of  which  were  prevented  from  rotating  but  were  free  to  move 
axially.  It  was  found  subsequently  that  better  agreement  between  experimental 
results  and  theoretical  predictions  was  obtained  when  the  dynamic  plastic 
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bending  moment  was  calculated  using  a  yield  stress  given  by  Manjoine  [8] 
corresponding  to  a  mean  value  of  strain-rate.  Symonds  and  Mentel  [9]  ex¬ 
amined  the  influence  of  axial  restraints  on  beams  loaded  impulsively  and 
predicted  final  deformations  which  were  considerably  smaller  than  those 
expected  from  the  corresponding  simple  beam  solution  when  deflections  of 
the  order  of  the  beam  depth,  or  larger,  were  permitted.  Humphreys  [10] 
conducted  some  experimental  tests  and  confirmed  this  prediction,  indica¬ 
ting  that  the  simple  theory  is  only  useful  for  maximum  deflections  of  mag¬ 
nitudes  up  to  the  order  of  the  beam  thickness.  In  reference  [11]  Florence 
and  Firth  report  the  results  of  some  experiments  in  which  pinned  and 
clamped  beams  without  axial  restraints  were  subjected  to  uniformly  distri¬ 
buted  impulses  large  enough  to  cause  considerable  plastic  deformation.  It 
was  found  that  a  rigid-plastic  analysis,  which  disregarded  strain-rate  ef¬ 
fects  entirely,  but  included  strain-hardening  in  an  approximate  manner  dur¬ 
ing  the  second  stage  of  motion,  gave  somewhat  better  agreement  with  the  ex¬ 
perimental  results  than  a  simple  rigid-plastic  analysis. 

Recently,  Nonaka  [12]  studied  the  behavior  of  clamped  beams,  with  re¬ 
straints  against  axial  displacements  at  the  supports,  when  an  attached  mass 
in  the  center  was  subjected  to  large  transverse  dynamic  loads.  It  was  ob¬ 
served  that,  except  when  the  attached  mass  was  small,  a  major  portion  of  the 
deformation  occurred  under  a  one  degree  of  freedom  mode  in  which  the  two 
halves  of  the  beam  rotate  about  the  supports.  Consequently,  the  mode  ap¬ 
proximation  method  of  Martin  and  Symonds  [13],  which  has  been  used  by 
Symonds  [14]  to  study  a  cantilever  beam  loaded  impulsively,  was  utilized  in 
order  to  estimate  the  effects  of  strain-rate  sensitivity,  elasticity,  and 
load  duration  on  the  final  deformation. 
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Apart  from  the  numerical  work  of  Witmer,  Balmer,  Leech  and  Plan  [15] 
the  combined  influence  of  strain-hardening  and  strain-rate  sensitivity  on 
the  permanent  deformations  of  rigid-plastic  beams  with  immovable  supports 
arising  from  uniform  impulses  has  not  been  studied  and  is  taken,  therefore, 
as  the  subject  of  interest  in  this  article.  In  order  to  retain  the  attrac¬ 
tive  simplicity  of  rigid-plastic  analyses ,  an  attempt  is  made  to  develop  a 
method  which  is  mathematically  simple  yet  at  the  same  time  sufficiently 
accurate  to  be  worthwhile  exploring  the  possibility  of  using  it  to  analyze 
the  behavior  of  more  complex  structures. 


2 .  General  Equations 

An  expression  for  the  axial  strain  at  any  position  in  a  beam  is  derived 
in  this  section  in  terms  of  the  membrane  or  axial  strain  and  curvature  of 
the  mid-surface.  The  equilibrium  equations  for  a  uniform  beam  loaded  dy¬ 
namically  are  then  derived  and  recast  into  a  form  convenient  for  later  use. 


2.1  Axial  Strain 


It  may  be  shown  that 


ds  =  ndx 


where 


n  =  {w'2  +  (l  +  u')2}1/2  , 


<  >'  >  • 


and  the  remaining  quantities  are  defined  in  Fig.  1. 
Thus , 

e  =  n  -  1 


(1) 

(2) 


(3) 
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is  the  axial  strain  of  the  mid-surface  of  the  beam  which  becomes 


s  =  u’  +  + 


when  n  is  expanded. 

The  radius  of  curvature  "R"  of  the  mid-surface  of  the  beam  is,  by 
definition. 


3*  * 


which,  using  the  geometrical  identity 


tan* 


from  Fig.  1(a)  and  equation  (1)  yields 


1  _  — w" ( 1  +  u* )  +  w'u" 


If  it  is  assumed  that  plane  cross-sections  remain  plane  during  deforma¬ 
tion  and  merely  rotate  about  the  mid-surface  of  the  beam,  then  it  may  be 
shown  that  the  axial  strain  at  distance  "z"  is 


where 


and  z  is  defined  in  Fig.  1(b). 


2.2  Equilibrium  Equations 


=  e  +  zk 


1  +  e 
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The  equations  of  equilibrium  for  the  element  of  the  beam  illustrated  in 


Fig.  1(c)  can  be  written  in  the  form 
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9x  ~  R 


yur|Cos<j>  + 


pw^sinij)  +  bn  =  0 


(8) 


9Q  .  Nn  . 

9x  R 


yiinsin<j)  +  yWncos<()  +  qn 
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(9) 


9M 

9x 


+ 


Qn 
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(10) 


provided  any  rotatory  inertia  effects  are  ignored  and 


<  >  =  lt  <  >  • 


In  order  to  simplify  the  theoretical  analyses  of  beams  and  strings 
loaded  dynamically,  the  displacement  u  and  acceleration  ii  are  usually 
disregarded  [9,  etc.].  Thus,  let 


(ID 


which  allows  equations  (2)  and  (4)  to  be  written 

o  1/2 

n  =  (l  +  w'2) 

and 


respectively. 

If  attention  is  restricted  henceforth  to  small  strains ,  then 


(12) 


(13) 


n  =  1  ,  sin4>  =  -w'  ,  cos<J>  =  1  , 


and  equations  (5)  and  (7)  give 


tc  =  -w" 


(14) 


Equations  (8)- (10)  may  now  be  rewritten 
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dm 

dx 


4w'n 

H 


(19) 


The  first  term  on  the  right  hand  side  of  equation  (19)  is  due  to  an  ex¬ 
ternal  load  k  ,  the  second  term  yw  is  an  inertia  force  while  the  last  term 
containing  w'  introduces  the  membrane  forces  arising  from  axial  restraints. 


3.  Yield  Condition 

It  is  indicated  in  Fig.  2  that  an  exact  yield  curve  relating  the  dimen¬ 
sionless  bending  moment  m  and  membrane  force  n  according  to  the  Tresca 
shear  stress  criterion  lies  everywhere  inside  a  square  having  sides  of  mag¬ 
nitude  2  ,  while  a  square  with  sides  of  length  1.236  lies  everywhere  in¬ 
side  the  exact  yield  curve.  In  order  to  estimate  the  accuracy  of  this  ap¬ 
proximate  linear  yield  condition,  which  can  be  selected  to  bound  the  exact 
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yield  condition,  as  illustrated  in  Fig.  2,  the  particular  case  of  a  rigid, 
perfectly  plastic  beam  loaded  impulsively  will  be  studied  and  the  results 
compared  with  the  final  deformations  predicted  by  Symonds  and  Mentel  [9], 
who  used  an  exact  yield  curve. 

4.  Rigid,  Perfectly  Plastic  Beam  Subjected  to  a  Uniformly  Distributed 
Impulse 

4.1  Simple  Bending  Solution 

The  equilibrium  equation  (19)  for  this  case  becomes 

=  !L_  »!»  (20) 

dx2  "o  n2 

if  the  beam  is  unloaded  (k  =  0)  and  there  are  no  axial  restraints  (w*  =  0) 
against  deformation. 

In  accord  with  experimental  evidence  and  previous  analyses ,  it  is  as¬ 
sumed  that  two  traveling  hinges,  each  of  radius  |p(t)|  ,  originate  from  the 
supports  at  t  =  0  and  travel  inwards  towards  the  center  of  the  beam  during 
a  first  phase,  while  throughout  a  second  stage,  they  remain  stationary  at 
x  =  0  until  all  the  initial  kinetic  energy  is  dissipated  as  plastic  work. 

First  Stage 

A  velocity  profile , 

w  =  V  for  0  S  x  $  p(t) 

o 

and 

V  (L  -  x) 

w  =  - z -  for  p(t)$x$L, 

L  -  p 


(21) 

(22) 
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where  Vq  is  the  initial  uniform  velocity  of  the  beam,  is  consistent  with 
the  normality  requirement  associated  with  the  yield  condition  m  =  1  . 

If  the  time  derivative  of  (22)  is  substituted  into  (20)  and  the  result¬ 
ing  expression  integrated  twice  with  respect  to  x  ,  one  obtains 

pV 

gg—  (L-p)p+l=0  ,  (23) 

o 

for  a  beam  which  is  simply  supported  at  its  ends.  Integrating  (23)  with  re¬ 
spect  to  time  yields 


t  =  —  (L2  -  2Lp  +  p2)  , 
L2 

where 

_  "V2 

T  "  6M 
o 


(24) 


(25) 


and  the  requirement  that  p  =  L  when  t  =  0  has  been  satisfied. 

The  maximum  displacement  at  the  center  of  the  beam  at  the  end  of  the 

first  stage  is  V  t  . 

o 


Second  Stage 

If  a  linear  velocity  profile  of  the  form 

*  =  W(t )  (-L-~  — )  (26) 

is  selected  for  the  second  stage,  then  it  may  be  shown,  using  (20)  and  (26) 
and  matching  the  displacements  and  velocities  at  t  =  t  with  those  at  the 
end  of  the  first  stage,  that 


where 


w  > 
m  _  1 

H  =  3 


(27) 


UV|L2 

M  H 
o 


X 


(28) 
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and  w^  is  the  permanent  displacement  at  the  center  of  the  beam. 

Equation  (27)  is  the  same  result  as  equation  (4)  of  reference  [9]. 


4.2  Solution  with  Axial  Restraints 


The  influence  of  membrane  forces,  which  arise  due  to  axial  restraints, 
must  now  be  included  in  the  equilibrium  equations.  Thus,  equation  (19)  be¬ 
comes 


dm 

dx 


o 


U  32w 

M 

o  3t2 


3W 

3x 


(29) 


First  Stage 

If  the  mechanism  of  deformation  for  this  case  is  assumed  to  be  similar 
to  that  for  the  previous  one ,  then  the  velocity  profile  is  given  by  equations 
(21)  and  (22),  and  when  0  $  x  $  p  . 


w  =  V  t 
o 


(30) 


while 


w  =  V  t(x)  + 
o 


t(x) 


3W  H+- 

st  dt  1 


if  n  £  x  <  L 


(31) 


Since  one  might  expect  bending  moments  to  dominate  over  the  action  of 
membrane  forces  during  the  first  stage,  then  the  form  of  the  time  function 
t(p)  for  wf  f  0  should  be  somewhat  similar  to  equation  (24)  which  was  ob¬ 
tained  by  disregarding  wf  in  (19).  Therefore,  it  is  assumed  that. 


t(p)  =  tQ(L2  -  2Lp  +  p2)  , 


(32) 


where  t  is  a  constant  to  be  determined  later, 
o 

It  may  be  shown  using  equations  (21,  22,  30,  31,  32)  that 


w'  =  0  ,  when  0  $  x  $  p 


w'=2Vt(p-x),  if  p  $  x  $  L  .  (34) 

o  o 

Equation  (29)  can  be  rewritten  with  the  aid  of  (22)  and  (34)  in  the  form, 

dm  _  yV°^  tr  x2  _  p2.  8Vo  /_  ^  /ocX 

ST  -  ~  ~  ■  ;  O.X  -  —  -  Le  ♦  — )  -  -g—  (P  -  »>  ,  (35) 

M  (L-p )z 
o 

provided  p£xSL,n  =  l  and  o  $  m  $  1  . 

Integrating  (35)  gives 

"V  ,Lx2  x3  T  p2x  .  Lp2  p3^  8Vo  x2  P2\  .  i 

M  „  n  \2  -  I"  -  ^  +  —  +  —  -  I" >  -  -TT-  <«---—)»  1  - 

MV  L-p ; 

(36) 

where  the  constant  of  integration  has  been  determined  from  the  requirement 
that  m  =  1  at  x  =  p  ,  which  in  turn  can  be  obtained  by  solving  equation 
(29)  using  (21)  and  (33). 

Now  for  a  beam  which  is  supported  at  pinned  ends ,  m  =  0  at  x  =  L  . 
Thus  (36)  gives 


pv  T  3  3 

- 2—  (i-  -  I,2p  *  Lp*  -  §-) 

2M  t  (L-p)3  d  6 


8V  t  T  2  2 

-P  (pL  -  f  -  fl)  +  1  =  0 


which,  when  p  =  0  ,  yields 


1  {.1+  (1  +  i|)1/2}  , 


where  6,  =  V  t,  is  the  maximum  deflection  (x  =  0)  at  the  end  of  the  first 
JL  O  1 

stage ,  and  from  ( 32 ) ,  t^  =  tQL2  ,  when  p  =  0  . 
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Second  Stage 

If  a  linear  velocity  profile  of  the  form  described  by  equation  (26)  is 

selected  for  the  second  stage,  then  it  is  straightforward  to  show,  when  using 

the  yield  condition  illustrated  in  Fig.  2,  that  the  final  displacement  of  the 
beam  at  x  =  0  is 

w  1/2  1/2 

-jp  ■  |  +  2X  ♦  J  (1  +  %>  )  -  1]  (37) 

The  final  deflections  at  the  center  of  rigid,  perfectly  plastic  beams 
loaded  with  uniformly  distributed  impulses  are  given  by  equations  (27)  and 
(37)  and  compared  in  Fig.  3  with  the  results  of  Symonds  and  Mentel  [9],  It 
is  evident  that  the  "upper"  and  "lower"  bounds  lie  between  the  two  predic¬ 
tions  of  Symonds  and  Mentel  [9]  with  and  without  a  string  phase,  while  the 
deflections  forecasted  by  the  simple  bending  solution  given  by  equation  (27) 
are  considerably  larger  even  for  very  small  values  of  A  .  Strictly  speak¬ 
ing,  the  curves  designated  "upper"  and  "lower"  bounds  in  Fig.  3  are  not 
upper  and  lower  bounds  in  the  accepted  sense  since  they  are  based  solely  on 
the  fact  that  the  deflections  designated  "upper"  bound  were  calculated  using 
a  yield  surface  which  lay  on  or  outside  the  exact  one ,  while  those  termed 
"lower"  bound  were  evaluated  using  a  yield  surface  which  lay  everywhere  on 
or  inside  the  exact  one. 

5 .  Constitutive  Equations 
5.1  Strain-Rate  Sensitivity 

It  is  well  known  that  the  initial  yield  stress  [8,  16,  etc.]  of  many 
materials  increases  with  increase  of  strain-rate  and  that,  furthermore,  it 
is  important  to  take  account  of  this  effect  when  analyzing  the  dynamic 
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behavior  of  cantilevers  and  beams  [4,17].  Cowper  and  Symonds  [18]  observed 
that  a  constitutive  equation  of  the  form 


1  ♦  <§> 


1/p 


(38) 


could  be  fitted  to  the  data  of  Hanjoine  [8]  provided  D  =  40.4  sec.  ^  and 
p  =  5  . 


On  account  of  its  attractive  simplicity,  the  constitutive  equation  (38) 
has  been  used  successfully  by  a  number  of  authors  to  solve  a  variety  of  prob¬ 
lems  [5,19,20]. 

If  an  element  of  the  beam  illustrated  in  Fig.  1  is  made  from  a  rigid, 
rate-sensitive  material  described  by  the  constitutive  equation  (38),  then 
it  may  be  shown  that 

2 


m  =  1  - 


(e  +  j  k) 


p+1 

P 


,H  • 

( J  K 


HD 


,1/P*  (P+1) 


+ 


a 2. 


H2ic2D1/p(2p+l) 


2p+l 

ic)  P  +  k  - 


_ !tei _ 

H2<2D1/P(p+1) 


Eli 


e) 


(39) 


provided  plane  cross-sections  remain  plane  and  during  deformation  merely 
rotate  about  the  mid-surface  of  the  beam,  and  that  the  directions  of  m  and 
n  correspond  to  those  of  M  and  N  in  Fig.  1. 
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When  m  =  0  and  k  =  0  ,  it  may  be  shown  that 

.  1/p 

n  =  1  +  (jj)  ,  with  £  i  0 


(40) 


while,  for  n  =  0  and  £=  0  , 


0  „.  1/p 

m  =  i +  2^?r  (#  »  with  *  *  ° 


(41) 


In  order  to  simplify  the  rather  complicated  constitutive  equation  (39), 
let  us  proceed  in  a  manner  somewhat  similar  to  that  described  previously  in 
section  3  for  the  rigid,  perfectly  plastic  case.  It  is  assumed  that  a  lin¬ 
earized  yield  curve  will  grow  as  illustrated  in  Fig.  4(a)  such  that  the  be¬ 
havior  suggested  by  (40)  occurs  along  A'B'  with  0  $  m  ^  1  ,  while  equa¬ 
tion  (41)  describes  the  behavior  of  C'D’  with  0  $  n  $  1  . 

Perrone  [20,21]  has  shown,  for  some  simple  structures  loaded  impulsively 
and  made  from  a  strain-rate  sensitive  material,  that  excellent  agreement  with 
exact  solutions  may  be  obtained  when  utilizing  a  strain-rate  insensitive  ma¬ 
terial  with  a  constant  yield  stress  equal  to  the  initial  dynamic  yield  stress. 
This  observation  permits  considerable  simplification  of  subsequent  analyses 
but  is  only  valid  for  large  values  of  p  (4  or  5)  when  most  of  the  kinetic 
energy  is  dissipated  before  the  stress-strain-rate  point  departs  appreciably 
from  its  initial  position. 

5.2  Strain  Hardening 

Cowper  and  Symonds  [18]  and  others  [11,  etc.]  have  suggested  a  linear 
strain-hardening  relation  of  the  form. 


(42) 
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where  "E"  is  the  Elastic  Modulus  and  "r”  may  be  interpreted  as  the  ratio 

of  the  slopes  of  the  elastic  and  plastic  portions  of  the  stress-strain  curve. 

* 

If  plane  cross-sections  are  assumed  to  remain  plane  and  merely  rotate 
about  the  mid-surface  of  the  beam  during  deformation,  then  it  may  be  shown, 
using  equation  (42),  that 

m  =  1  -  (n  -  ve)2  +  — -  ,  (43) 

provided  "m"  and  "n"  act  in  the  directions  of  "M"  and  "N"  indicated 

in  Fig.  1,  and  v  =  —  . 

°  ’  a  r 

o 

When  m  =  k  =  0  ,  it  may  be  shown  that 

n  =  1  +  ve  ,  provided  e  £  0  ,  (44) 


while,  if  n  =  e  =  0  ,  and  k  z  0  , 


m 


1  + 


vHk 

3 


(45) 


A  procedure  similar  to  that  outlined  previously  in  sections  3  and  5.1  for  the 
rigid,  perfectly  plastic  and  rigid,  strain-rate  sensitive  cases  will  be  adop¬ 
ted  here  in  order  to  simplify  the  constitutive  equation  (43)  for  a  rigid, 
strain-hardening  material.  It  is  now  assumed,  therefore,  that  a  linearized 
yield  curve  illustrated  in  Fig.  4(b)  will  grow  as  indicated  so  that  the  side 
A'B'  is  described  by  equation  (44)  with  0  $  m  $  1  ,  while  (45),  with 
0  <  n  $  1  ,  defines  the  behavior  of  side  C'D*  . 


5 . 3  Combined  Strain-Rate  Sensitivity  and  Strain  Hardening 
Symonds  [14]  and  Perrone  [22]  suggested  that 


2-  =  f(e)  g(e) 
a 


(46) 
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could  be  used  to  analyze  structures  loaded  dynamically,  where  f(e)  and 
g(e)  a,re  strain-rate  sensitivity  and  strain-hardening  relations,  respec¬ 
tively.  It  is  well  known  that  strain-hardening  of  some  materials  decreases 
with  increase  in  strain-rate,  and  strictly  speaking  a  stress-strain-strain 
rate  law  cannot,  therefore,  be  written  in  the  product  form  of  equation  (46) 
with  f(e)  and  g(e)  uncoupled.  However,  in  order  to  retain  mathematical 
simplicity,  it  is  assumed  that  the  combined  effect  of  strain-hardening  and 
strain-rate  sensitivity  could  be  considered  in  the  manner  suggested  by  (46), 
which  using  (38)  and  (42)  may  be  rewritten, 

~  =  {i  +  (§)  j(l  +  ve)  (47) 

Thus,  from  equations  (40,  41,  44,  45,  47)  it  is  evident  that,  when 
e  £  0  ,  t  i  0  ,  <  =  0,  ic  =  0  and  0  $  m  S  1  ,  then 

.  1/p  .  1/p 

n  =  1  +  ve  +  (^-)  +  ve  (^-)  (48) 

and  when  e=0,  £  =  0  ,  £^0,  k^O,  and  0  $  n  $  1  ,  then 


m 


1 


+ 


2p+l 


(Hi)1/P 

'2D' 


,  vHk  .  2vHicp 
+  3  +  3(2p+l)  (2D* 


(49) 


6 .  Influence  of  Strain-Hardening  and  Strain-Rate  Sensitivity  on  Impulsively 
Loaded  Beams 

In  order  to  examine  the  influence  of  strain-hardening,  strain-rate  sen¬ 
sitivity,  and  the  combined  effect  of  both  on  the  large  final  deformations  of 
beams,  the  constitutive  equations  developed  in  the  previous  section  will  be 
used  to  analyze  a  rigid-plastic  beam  with  axial  restraints  which  is  simply 
supported  at  pinned  ends  and  loaded  with  a  uniformly  distributed  impulse. 
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6.1  Strain-Rate  Sensitivity 
First  Stage 


Equations  (13)  and  (14)  give 


e  =  w'w'  ,  and  k  =  -w"  .  (50) 


Thus,  it  may  be  shown  when  employing  the  mechanism  of  deformation  for 
the  rigid,  perfectly  plastic  beam  analyzed  in  section  4.2  that  the  velocity 
profile  described  by  equations  (21)  and  (22)  is  consistent  with  the  yield 
condition  (40),  provided  the  implication  of  a  discontinuity  in  curvature  at 
the  traveling  hinge  is  disregarded.* 

Using  the  velocity  profile  described  by  (21),  the  strain-rate  sensitive 
relation  (40),  and  equation  (33)  in  order  to  solve  the  equilibrium  equation 
(29),  yields 


m  =  1  ,  for  0  $  x  $  p 


When  p  $  x  $  L  ,  the  corresponding  equations  give 


^  it  is  assumed  that  the  discontinuity  of  curvature  given  by  the  differ¬ 
ence  of  the  appropriate  derivatives  of  equations  (21)  and  (22)  can  be  re¬ 
placed  by  a  continuous  change  of  curvature  across  an  annulus  of  width  2H 
then  it  may  be  shown  that 

P  =  -r—  ,  approximately 

where.  "P"  is  the  ratio  of  energy  dissipated  at  a  traveling  hinge  during 
the  first  stage  to  the  corresponding  loss  of  kinetic  energy,  and  "a"  in¬ 
dicates  the  factor  by  which  "m"  is  increased  according  to  the  strain-rate 
sensitivity  relation  given  by  equation  (41). 

In  order  to  simplify  the  analysis,  average  values  have  been  used  for  the 
speed  of. the  traveling  hinge  and  the  radius  of  curvature  across  the  annulus. 
P  ,  as  might  be  expected  has  rather  large  values  when  X  is  small  but  for  a 
mild  steel  beam  with  L/H  =  12  (L  =  12  in.)  and  X  =  100  ,  P  =  0.06 
while  for  X  =  800  ,  P  =  0.01  .  ’ 
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-  mV 


,  Lx2  x3 


(  U.A.  A.  T  P  2X  L  p  2  P  3\ 

m  =  -  (-= - - - Lpx  +  +  -= - «— )  + 

2M  t  (L-P)3  26  223 


o  o 


2N  V  t  f 

O  O  O  3 

— m - : 

O 

m 


2  2 

D 

— — —  —  Dv  4"  -  +  -  £■ 

2  M  2  2p+l 


\  D(L-p)y 


1/p 

(x  -  p)  v  +  1 


) 


(51) 


where  the  condition  m  =  1  and  the  continuity  requirement  [m]  +  (5  [m  *  ]  =  0 
have  been  satisfied  at  x  =  p  . 

If  the  beam  is  simply  supported  at  pinned  ends,  then  the  requirement 
that  m  =  0  at  x  =  L  allows  equation  (51)  to  be  rewritten  in  the  form 


2p+l 

Sg? t*  p 


(52) 


when  p  =  0  ,  and  where 


and 


(53) 


(54) 


Second  Stage 

Humphreys  [10]  and  Florence  and  Firth  [11]  have  conducted  experiments 
on  beams  loaded  impulsively,  from  which  it  is  evident  that  the  deflections 
occurring  at  the  end  of  the  first  stage  are  comparable  with  the  beam  depth. 

It  is  shown  in  references  [9,23,24]  that  membrane  forces  dominate  over  bending 
moments  when  deflections  are  of  the  order  of  the  thickness  of  rigid-plastic 
beams  and  plates  loaded  dynamically.  Furthermore,  the  mechanism  of  deforma¬ 
tion  suggested  by  equations  (21)  and  (22)  implies  that  two-thirds  of  the 
initial  kinetic  energy  is  dissipated  during  the  first  stage,  leaving  one- 
third  to  be  dissipated  in  the  final  stage. 
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In  view  of  the  foregoing  comments,  it  seems  reasonable  to  consider  that 
the  beam  behaves  as  a  membrane  during  the  second  stage,  the  equilibrium  equa¬ 
tion  of  which  may  be  obtained  from  (16)  with  Q  =  0  , 


where  from  (40) 


w"  =  -U-  w 
N  n 
o 


n  =  1  +  ( 


and  m  =  0  . 


(55) 


(56) 


It  is  shown  later,  for  the  rigid,  perfectly  plastic  case,  that  equation 
(55)  predicts  final  results  almost  identical  to  those  of  equation  (37).  In  order  to 
simplify  the  solution  of  (55)  use  will  be  made  of  Perrone’s  [20]  observations 
noted  previously  in  section  5.1.  Thus,  at  the  end  of  the  first  stage  when 
p  =  0  ,  equations  (22,  34,  50)  give 


(57) 


which,  when  substituted  into  (56),  yields 


/2V  2 t  *\1/P 


the  average  value  of  which  is 


(59) 


It  is  more  convenient  to  use  "y"  measured  from  the  left-hand  support 
so  that  y  =  0  ,  2L  at  the  supports  and  y  =  L  at  the  beam  center;  and  in 
order  to  remove  the  non-linearity  arising  in  (55)  due  to  (58),  a  solution 
will  be  sought  using  nD  instead  of  n  . 
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Now  the  beam  finally  comes  to  rest  at  t  =  ,  where  t^  =  t^  +  T^  . 
Thus,  equating  the  time  derivative  of  (60)  to  zero  and  using  (64)  and  (65) 
gives 


kirc-T.j. 

/  X\  X  V 

c°s  (  — ~ 


±2cRLto 

(1  +  4c2L2t  2)1/2 

K  O 


where  +  sign  is  to  be  used  when  k  =  1,5,9  and  -  sign  when  k  =  3,7,11,... 

Making  use  of  equations  (60,  64,  65,  66,  67)  the  deflection  of  the  beam 
at  T  =  Tj.  and  y  =  L  can  be  expressed  in  the  following  form. 


where 


and  from  (59), 


»  -I  K  ,  J-/  * 

TT  =  (!m7  *  U42)  * 

IT3  R 


*  =  l  i_ 

k=l,3,5  k3 


(68) 


(69) 


(70) 


Equation  (68)  with  a  =  0  reduces  to  the  rigid,  perfectly  plastic  case 
and  predicts  final  deflections  at  the  center  of  the  beam  which  when  plotted  in  Fig. 
3  are  almost  coincident  with  the  results  of  equation  (37)  for  the  "upper"  bound 
case.  Thus  waiving  the  requirement  that  m  should  be  continuous  with  re¬ 
spect  to  time  between  the  two  stages  of  deformation  leads  to  considerable 
simplification  with  no  concomitant  loss  of  accuracy  provided  that  the  L/H 
ratio  of  the  beam  is  not  too  small. 
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6.2  Strain-Hardening 
First  Stage 

It  may  be  shown  that  following  a  procedure  somewhat  similar  to  the  one 
outlined  in  section  6.1  for  the  strain-rate  sensitive  case,  but  using  (44) 
instead  of  (40),  gives 

4YA4  +  4A2  +  A  -  ^  =  0  ,  (71) 

b 

where 

Y  =  v(£)  (72) 

and  the  implications  of  a  curvature  discontinuity  at  the  traveling  hinge 
have  been  disregarded.  In  fact,  an  analysis  similar  to  the  one  described 
in  the  footnote  of  section  6.1  indicates  that  P  is  of  the  same  order  of 
magnitude  as  in  the  rate-sensitive  case. 


Second  Stage 


It  is  assumed  that  the  beam  behaves  like  a  string  throughout  the  second 
stage  since  this  simplification  for  the  rigid,  perfectly  plastic  case  has 
been  shown  previously  to  lead  to  final  deflections  which  are  almost  the  same  as 
those  predicted  by  equation  (37). 

At  the  end  of  the  first  stage,  t  =  t.  ,  and  from  equations  (34)  and  (44) 

n  =  1  +  2vV2t2x2  ,  (73) 

o  o  ’  v  ' 


the  mean  value  of  which  is 


n 


H 


=  1  + 


2vV2t2L2 
o  o 


3 


(74) 
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If  the  equilibrium  equation  (55)  was  solved  using  n  =  n^  ,  then  any 
strain-hardening  occurring  during  the  second  stage  would  be  disregarded. 
Thus,  in  order  to  consider  this  additional  strain -hardening  in  an  approxi¬ 
mate  manner,  the  final  strains  corresponding  to  n  =  nH  will  be  found  by 
solving  equation  (55).  This  then  allows  a  more  realistic  estimate  of  n 
to  be  made  by  using  the  average  of  the  value  of  n  corresponding  to  this 
final  strain  and  n^  . 

Matching  the  displacement  and  velocity  of  the  beam  at  the  end  of  the 
first  stage  to  those  at  the  beginning  of  the  second  in  the  manner  described 
previously  in  section  6.1  gives 


16V  L  1/2 

- —  (1  +  4c2L2t2) 


,(*n.) 

v  21/ 


,(*!£) 

V  21/ 


k=l,5 ,9 


k=3,7,ll 


- —  (1  +  4c2L2t2) 

«2ch  H  ° 


I 

k=l  ,5  ,9 


><Ht> 


k=3 ,7 ,11 
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It  is  evident  from  (44)  that  a  membrane  force 


np  =  1  +  ve 


(79) 


corresponds  to  the  final  strain  e  given  by  equation  (77)  which  was  derived 
assuming  no  strain-hardening  during  the  second  stage. 

Now,  assuming  that  any  strain-hardening  during  the  second  stage  may  be  ac¬ 
counted  for  by  taking  the  average  of  membrane  forces  n^  at  T  =  0  (t  =  t^)  and 
np  at  T  =  Tp  (t  =  t^  +  Tp)  ,  then 


nu  ,  8vV  1  . 

n .  =  -y  +  \  +  - —  (—  +  4L2t2)  a  , 

A  2  2  ,4  „  2  o 


(80) 


7T  C 


H 


which  gives,  finally 


(  X  +  4a2) 


1/2 


3  4nA 
ira  A 


(81) 


6.3  Combined  Strain-Rate  Sensitivity  and  Strain-Hardening 
First  Stage 

If  procedures  which  were  developed  for  the  analyses  of  the  first  stages 
in  sections  6.1  and  6.2  are  followed,  but  using  equation  (48)  instead  of  (40) 
or  (44),  then  it  may  be  shown  that 

4p+l  2p+l 

&  *  p  -jSr  4  p  +4Yi-t4s2  +  i.x  =  0  (82) 


from  which  the  rigid-plastic  result  and  equations  (52,71)  can  be  obtained  by 
putting  a  =  y  =  0  ,  y  =  0  and  a  =  0  ,  respectively. 
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Second  Stage 

If  it  is  assumed  that  membrane  forces  alone  are  important  during  the 
second  phase  of  deformation  and  that  strain-rate  sensitivity  and  strain¬ 
hardening  may  be  accommodated  in  the  manner  described  in  section  5.3,  then 


where 


and 


t  '  7  <55- +  “*2>  t 

TJ3  C 


n  =  (i  +  ES*_)(i  +  2|i+  Mi.  + 

C  ptl  3  ir4n  ,* 

o 


n  =  (1  +  22^_)(1  ♦  2lfi) 
o  p+1  3 


(83) 


(84) 


(85) 


Equation  (83)  reduces  to  the  rigid,  perfectly  plastic  case  and  equations 
(68)  and  (81)  when  a  =  y  =  0  ,  y  =  0  and  a  =  0  ,  respectively. 


Discussion 

It  is  clearly  evident  from  Fig.  5,  which  is  plotted  using  values  from 
equation  (68)  for  the  strain-rate  sensitive  case,  that  it  is  important  to  in¬ 
clude  strain-rate  effects  when  estimating  the  permanent  deformations  of  beams 
loaded  impulsively.  It  is  interesting  to  note  for  given  values  of  the  im¬ 
pulse  parameter  X  ,  L/H  ratio,  and  material  properties,  that  the  maximum 
deflection-to-thickness  ratio  is  smaller  for  smaller  beams.  This  "size  effect," 
which  has  been  noted  previously  by  Symonds  [14]  and  others  and  follows  from  the 
form  of  equations  (52)  and  (53),  is  indicated  in  the  Appendix. 
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Fig.  6,  which  illustrates  the  influence  of  strain-hardening  given 
by  equation  (81),  shows  that  strain-hardening  is  particularly  important 
for  beams  having  small  L/H  ratios. 

The  results  for  situations  where  strain-hardening  and  strain-rate 
effects  exist  simultaneously  are  given  by  equation  (83)  and  plotted  in 
Figs.  7-9.  These  curves  indicate  that  including  either  strain-hardening 
alone  for  beams  with  small  L/H  ratios ,  or  strain-rate  sensitivity  alone 
for  physically  small  beams ,  or  either  for  other  beams ,  gives  results 
which  compare  quite  favorably  with  the  combined  case.  It  also  appears 
from  Figs.  8  and  9  which  are  plotted  for  steel  and  aluminum  beams  having 
the  material  properties  listed  in  Table  I  and  other  calculations  that 
strain-rate  sensitivity  and  strain-hardening  are  not  important  for 
physically  large  beams  having  large  L/H  ratios . 

Conclusions 

A  simple  method  is  herein  suggested  for  estimating  the  combined  influ¬ 
ence  of  strain-hardening  and  strain-rate  sensitivity  on  the  permanent  defor¬ 
mation  of  structures  loaded  dynamically.  In  order  to  assess  the  predictions 
of  the  linearized  rigid,  strain-rate  sensitive,  strain-hardening  constitutive 
equation  (48),  a  study  has  been  made  of  the  behavior  of  a  beam  which  is 
loaded  impulsively  and  supported  at  its  ends  by  immovable  frictionless 
pins.  It  is  evident  from  Figs.  5-9  that  when  considering  strain-hardening 
alone  for  beams  with  small  L/H  ratios,  or  strain-rate  sensitivity  alone 
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for  physically  small  beams,  or  either  for  medium  ones,  then  permanent 
deflections  are  predicted  which  compare  rather  favorably  with  those 
given  for  the  same  value  of  \  by  an  analysis  retaining  their  combined 
influence.  Moreover,  the  results  suggest  that  it  is  not  necessary  to 
include  either  strain-hardening  or  strain-rate  sensitivity  for  physically 
large  beams  having  large  L/H  ratios. 

It  is  not  possible  to  compare  the  theoretical  predictions  of  this 
article  with  experimental  results  since,  to  the  author's  knowledge,  no 
test  data  have  been  published  for  the  particular  case  analyzed.  However, 
it  is  rather  encouraging  to  note  that  Humphreys  [10]  in  a  study  of 
clamped  mild  steel  beams  loaded  impulsively  recorded  permanent  deflec¬ 
tions  which  had  the  same  proportion  of  the  values  predicted  for  a  rigid, 
perfectly  plastic  beam  as  forecasted  here  for  the  corresponding  pinned 
case. 

It  is  thought  that  the  method  suggested  herein  could  in  principle, 
at  least,  be  used  to  analyze  beams  having  other  support  and  loading 
conditions  and  extended  in  order  to  examine  the  behavior  of  plates  and 
shells,  though  it  is  felt  that  some  supporting  experimental  results  are 
required  in  order  to  assess  the  validity,  or  otherwise,  of  the  various 
approximations  made  in  the  theory. 
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Appendix 


Size  Effect 


X  =  constant 


Consider  two  beams  each  of  unit  width  which  have  the  same  material 
properties  P  ,  D  ,  p  ,  and  L/H  ratio  and  X  but  different  Vq  , 


L  ,  and  H  . 


Now,  from  equation  (28)  X  * 


2  2  2  2 
4yV  L  4pV  L 

o  _  o 

a  H3  o  H2 

o  o 


(i) 


Thus  if  X^  =  X2  and  L^/H^  =  ,  then 


V  ,  =  V  _ 
ol  o2 


(ii) 


Also 


o  = 


W  oH 
DL2 


\1/P 


hence , 


('al\P  VolL2 

Vi)  Vo2^1  ’ 


(iii) 


which  using  (ii)  becomes 


H, 


(iv) 


i.e.,  if  a2  >  al  ’  t*ien  H2  <  ^1  *  **  ev^^ent  from  Figs.  5  and  7-9 

that  smaller  deflection-to-thickness  ratios  are  obtained  as  a  increases 
for  a  given  value  of  X  .  Thus,  physically  smaller  beams  are  more  sensi¬ 
tive  to  strain-rate  sensitivity  than  larger  ones. 


Material 


Mild  Steel 


Table  I 


aQ,  psi  y/H  lb  sec2/!!!**  v  D,  sec  1  p 


30  000  0.000732  6  40.4  5 


Aluminum  6061-T6 


40  000 


0.000253 


2.85 


6500 
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Figure  2 
F igure  3 

Figure  4 

Figure  5 

Figure  6 
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Figure  9 


(a)  Deformed  Shape  of  Mid-Plane. 

(b)  Curvature  of  Beam. 

(c)  Forces  and  Moments  Acting  on  an  Element  of  the  Beam. 

-  Yield  Conditions  for  a  Rigid,  Perfectly  Plastic  Material. 

Impulsive  Loading  of  a  Simply  Supported  Rigid-Plastic  Beam 
with  Axial  Restraints. 

(a)  Rigid,  Strain-Rate  Sensitive  Yield  Condition. 

(b)  Rigid,  Strain -Hardening  Yield  Condition. 

Size  Effects  for  Constant  X  in  Impulsively  Loaded  Beams 
which  are  Made  from  a  Strain-Rate  Sensitive  Material  and 
Restrained  Axially. 

Influence  of  Strain -Hardening  Alone. 

-  Combined  Influence  of  Strain-Hardening  and  Strain-Rate 
Sensitivity. 

-  Combined  Influence  of  Strain-Hardening  and  Strain-Rate 
Sensitivity  for  Mild  Steel  Beams  with  the  Material  Charac¬ 
teristics  Listed  in  Table  I. 

-  Combined  Influence  of  Strain -Hardening  and  Strain-Rate 
Sensitivity  for  Aluminum  6061-T6  Beams  with  the  Material 
Characteristics  Listed  in  Table  I. 
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THEORY 
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